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Motivated by link-transformations of lattice gauge theory, a method for generating local unitary 
invariants, especially for a system of qubits, has been pointed out in an earlier work [M. S. Williamson 
et. at, Phys. Rev. A 83, 062308 (2011)]. The current work first points the equivalence of the so 
constructed transformations to the combined operations of partial transpose and realignment. This 
allows for a natural generalization to any system, with subsystems of arbitrary dimensions, in a 
straightforward manner and is shown to indeed generate local unitary invariants. Some properties 
of the resulting operators and simple consequences for tripartite higher dimensional states are briefly 
£f~) • discussed. 
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Entanglement, the remarkable nonlocal feature of quantum mechanics, has been extensively studied in the recent 
past 0. Owing mainly due to its role in quantum information protocols, for example in teleportation [2j, dense coding 
Q, and channel discrimination Q, as well as due to its presence in quantum algorithms that have a speedup over 
classical ones §. It also plays an increasingly prominent role in condensed matter physics @, and presumably also 
affects macroscopic observables such as magnetic properties in some solids @- 
i—,' However, to detect and measure entanglement in a general state, represented by a density matrix, is a difficult 
problem. It has been solved for the case of two qubit states (§l-[loj. bipartite states of a qubit and a qutrit [lfj, and 
C^' all bipartite pure states. There are various measures of entanglement, for example the von Neumann entropy of any 
> \ one subsystem in a bipartite pure state is such a measure, the concurrence introduced in measures entanglement 
between two qubits in a pure or mixed state, while negativity and log-negativity [Til fl2l| is invoked for a g eneral 
bipartite mixed state which uses the positive, but not completely positive, map of partial transposition [13| |. The 
^-i partial transpose was introduced to detect entanglement, and provides a necessary but not sufficient condition. One 
i— i' other such criterion uses the operation of "realignment" to be expanded on further below. 

Another approach to study entanglement, especially of a multipartite kind, is by studying the local unitary (LU) 
V/ invariants of the system which consists of invariants under arbitrary unitary transforms restricted to the individual 
£^ | subsystems. In as much as entanglement quantifies non-local properties, entanglement measures remain invariant 
under LU operations and hence the importance of their study. The spectra of the density matrix itself and the various 
reduced density matrices got by tracing out subsystems are such LU invariants. However, it helps to have invariants 
. that are polynomials in the entries of the density matrix [14| . and those whose physical interpretation in terms of 
entanglement is available. These invariants uniquely determine the orbit of the state under these local operations. 
For example three qubit pure states have five independent LU invariants (l5l HH ] excluding the trace. In a recent 
paper Williamson et.al [17| . inspired by lattice gauge theory, have given a method of generating these invariants by 
associating them with a closed path joining some or all the qubits where two consecutive qubits on the path are 
connected by a "link transformation" . More recently Szilard Szalay in [l8j has given index-free formulae for invariants 
• »"H | of k qudit system up to degree six by using graph-theoretic methods. 

In the present work, it is first shown that the link transformation in [l7j is unitarily equivalent to the combined 
operations of partial transpose [13] and realignment [l9[ (refer Eq. ([7])). Thus it is interesting that two rather 
independent operations on which entanglement criteria are based come together in the construction of local invariants. 
This result immediately suggests a way to generalize to a system of an arbitrary set of qudits, each not necessarily 
of identical dimensions. Such a generalization is then shown to be LU invariant as well. One of the advantages of 
this method is that it does not need the generalization of Pauli matrices in higher dimensions to get the invariants as 
required for the link transformation approach of fl7| . 

For the present, it helps to first summarize the central result of [171 ]. Consider a lattice of points a\, 02, . . . , 
where each point represents a qubit and N is the total number of qubits. The invariants are generated by various 
closed path (containing some or all the qubits) and two consecutive qubits are connected by a link transformation. For 
example consider a closed path as shown in Fig. (1) and qubits a\ and are lying on it then the link transformation 
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connecting them is denoted by S(a,2,ai) whose elements are given by 

S(a 2 ,ai) nm = ^tr[pi2 a% <g> < 2 ]. (1) 

Here p±2 is the reduced density matrix of qubits a\ and 0,2, o~i (I — 0, 1,2,3) are 2x2 matrices such that / = 
corresponds to the identity matrix and I = 1,2,3 corresponds to the standard Pauli matrices o~ x ,a y ,o~z respectively. 
In a similar way consecutive qubits are connected on the path by corresponding transformation matrices. The 
interpretation of £(02,01) was provided as m.2 — S(d2> fli)mi, where mi represented an arbitrary measurement 
operator in the space of subsystem 01, while 1112 was the resultant state of subsystem 02, all of these being written in 
the {ai} basis. 

Then for a given closed path, for example consisting of K qubits, not necessarily all distinct: {a}i< = 
(ax, d2, . . . , clk) , the quantity 

tr[5(ai, a K ) ■ ■ ■ S(a 3l a 2 )S(a 2 ,a 1 )] (2) 

was shown to be a LU invariant. Note that this equation is to be read from right to left. Different path coordinates 
ai can refer to the same subsystem, depending on the actual path. One can obtain invariants by taking various closed 
paths. In |17| a maximal set of five independent invariants of three qubit pure states were calculated in this manner. 

It is now shown that realignment of a partial transpose is a unitary transform of the link transformation S(a2, on). 
In the partial transpose (PT) operation, transposition is done only on one subsystem. As mentioned already, PT is 
a positive but not a completely positive map and can hence be used to detect entanglement. Consider a bipartite 
density matrix pi2, which refers to the reduced density matrix of the subsystems at points a\ and a2 of the path, and 
any orthonormal product basis in the corresponding state space, then the transposition on only the subsystem 

labeled 02, denoted as p^, is given by the matrix elements: 

(iMp%\j)\a) = (i\(a\ Pl 2\j)W). (3) 

The PT criterion states that if p\\ is negative then the state is entangled but otherwise it may or may not be 
separable, with separability guaranteed only for 2x2 and 2x3 systems This implies that positivity under PT 
is a necessary but not sufficient condition for separability. 

The other operation of interest is realignment, to which is associated what is called in literature as the computable 
cross norm (CCN) criterion [l9l. [20|. or simply as the realignment criterion. The corresponding operation on the 
density matrix p 12 , denoted as 1Z (P12) is given by: 

(i\{j\(K(p 1 2))\a)\P) = {i\{a\p 1 2\m. (4) 



The criterion is that if the state p\2 is separable then ||7£(/9i2)||i < 1, where ||M||i is the trace norm equal to trV MM^ 
[2lj . This condition is found to detect some bound entangled states, these being positive under PT and hence not 
being detected by the corresponding criterion [2(J H3| • 

Define the map P12 i-> 7?.(/5^), where both PT and realignment is effected serially. As both of these operations are 
merely permutations of the matrix elements, so is this combination. It is now shown that for two qubits, 5(01,02) is 
unitarily equivalent to IZlp^). Starting from Eq. ((T|) and using an arbitrary orthogonal product basis |i/3) (where i, 
/3 = 0, 1) and Einstein's summation convention (operative throughout the paper) one obtains the following: 

S{a 2l a 1 ) nm = -<*jS|pia <C ® 



(it3\p 12 \ja)(ja\o-%®o-%>\0) (5) 
(Pa\KO$i)\ji)<j\<r%\i)(a\o?\0). 



Define a 4 x 4 matrix U with elements: 



(ji\U\m) = -^0Vm|i>, (6) 

where i, j=0,l and ra — 0, 1, 2, 3. The matrix U is independent of the qubit label, and hence these have been omitted 
from the Pauli matrix symbols. Thus using Eq. ^ and the last expression of Eq. (JS|) gives the following equation: 

S(a 2 ,a 1 ) nm = (aP\U\n)(l3a\1l{p%)\ji)(ji\U\m). 
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On taking the complex conjugate of Eq. ([6]) and using the definition of U matrix, it can be shown that (a.p\U\n) — 
(n\W\f3a). Thus one obtains the advertised relation between the link transformation (S(a 2 ,ai)), and the operations 
of partial transpose and realignment as: 

S(a2, ai )=tfK(j$i)U. (7) 
The matrix U written explicitly in the standard basis is: 

/ 1 1 \ 

V i o o-ij 

which is the identity and the three Pauli matrices deformed into 4- vectors and arranged along columns in a matrix. It 
is readily seen that the matrix U is indeed an unitary matrix with the additional property that UU T = S, the swap 
operator. An invariant corresponding to the closed loop (a±, a 2 , . . . , clk) is then 

tv[S( ai ,a K ) ■ ■ ■ S(a 3 , a 2 )S(a 2 , ai )} = tr[U*1l(fi*%)U ■ ■ ■ U<<K(p%)UUtK(p%)U] 

Hence the link transformation matrices such as S(a 2 , a±) can be replaced by TZ(p 2 l), and referred to as transformation 
matrices themselves. That these are link transformations in their own right and are in fact amenable to the identical 
interpretations as S(a 2 , ai) in [ljj follows from expressing the measurement operator at subsystem a\ in the standard 
basis ejj, which are matrices with 1 at positions ij and elsewhere and i,j G {1,2}. That is the following is the 
equivalent to Eq. (jTJ) : 



(fc(p£)) 



= tr 



Pi 2 e^0(e-) T 1. (10) 



Generalizing beyond qubits to subsystems of higher dimensions di (qudits) it is shown in Proposition ([1} that the 
last term of Eq. (JSJ is a LU invariant, thus obviating the need to specially generalize link transformations such as 
S(a 2 ,ai). The quantity lZ(p 2 l) is in general a rectangular array of dimension d 2 x d\. The dimensions of the various 
reactangular arrays dovetail such that the final array on completion of the loop based at 1 is a square matrix of 
dimension d\. 

Proposition 1. Under local unitary operations, if {o}k = (ai, a 2 , . . . , ax) is a closed path, 

V({a} K ) ee H(j%*) ■ ■■Tl{pll)'R{p T 2 l) = (I7i ® U^TZ(p^) ■ ■ ■n(p? 2 )K(p 2 r [)(U 1 ® U?). (11) 

Proof. This follows immediately from the following observations: If 

P2i = (U 2 ®h)p 21 {Ul®h) (12a) 
and p 32 = (h®U 2 )p 2Z {h®ul), (12b) 

then n{pll) = {U 2 ®U* 2 yn{g{), (13a) 

K{pll) = n(p T £)(U 2 ®UZ), (13b) 

K{p2i) = (U 2 <Z>U^K(p 21 ), (13c) 

and K{ P32 ) = U{p 32 )(U;®U 2 ). (13d) 

Consider the case of only a local unitary U 2 acting on the second subsystem. Under this operation p\ 2 <— > p 2 \ and 
hence p 2 \ = (U 2 <E) h)p2i(U 2 It follows that: 

{P2l)ia;jfs = (U 2 t ) l ^(p 21 ) t > a - J > fs{U 2 )j' ,j ■ (14) 

Using the definitions of the realignment and partial transpose operations the above is rewritten as: 

(npzi)) ..= (uh,' {npii)) t (u 2 ) rj 

= (Ul)iA'(Ui)l j ,(K(pT{)). t . 

J ' J V /t'l'-Ba 
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CLj{ CL\ 0,2 O.3 



n(p T 2 \) = US{a 2 ,a 1 )U^ 

FIG. 1: (Color online) A closed path connecting some qudits of a total of K qudits is shown. Also shown is a fully retracing 
path on four qudits. Any two consecutive qudits on a path are connected by a transformation matrix, for example, qudits a\ 
and (22 are connected by 72-(pif)- The equality shown corresponds to the case of qubits. 



which is the the same as Eq. (|13aj) . Similarly, using Eq. (I12b[) leads to Eq. (I13b|) . While Eqs. (|13cll3d[) . derived in a 
similar way, for the case when there is only realignment and no partial transpose look the same, note the different 
ordering of the unitary matrices in them. Thus if Ui, U2, ■ ■ ■ Uk are local unitary operators acting on subsystems 1, 
2, . . ., K then their combined action on the successive bipartite states pi+u leads to the claim in Eq. ([TT]) . □ 

Thus it follows that the eigenvalues of V({cl}k) are LU invariants, as well as its trace. This provides the generaliza- 
tion of the link transformation approach to generate local unitary invariants by using realignment and PT operations 
i.e. replacing S(a,2, ai) by 72.(p^J), see Fig. ([1]). In fact the above proof does not need special properties of the Pauli 
matrices used in [17| . and obviates completely the need for their generalizations and hence presents a much simpler 
algorithm. 

The need for the joint operations of realignment and PT is interesting and if say one uses only realignment, it is 
easy to see that this does not in general lead to LU invariants. This follows from the right hand sides of Eqs. (|13cll3dp . 
where U2 <8> and 1/% <8> U2 appears. Thus if orthogonal transforms are used instead of the more general unitary 
ones, the operation of PT can be dispensed with. Alternatively if one restricts oneself to a real state space (such as 
the so-called rebit space [23j]) these will in fact be local invariants. 

Given a path {a}x and the associated operator V({ci}k), there are many different ways of associating LU invariants, 
even if one were to restrict to polynomials of the density matrix elements. For instance the coefficients of the 
characteristic polynomials of the operator are such invariants and these include the trace and determinant. That such 
coefficients are real follows from a property that is stated as part of Proposition ([2]) . Also if the path is fully retracing 
(there exists a number L such that the path is {01, . . . , at, Ol-i • • • ffli}, the retracing path shown in Fig. ([T]) has 
L = 4.) then V({o}k) is a positive operator itself. A necessary ingredient in the properties of the the operator is the 
swap operation. Consider the swap operator Si which acts on the product space of two identical dimensional spaces: 
®Hdx- Its action is given by (lm\Si\ij) — SijS m i. The swap operator is symmetric permutation matrix such that 
Si = I which implies further that it is unitary. Swap operators for other subsystems, Si are similarly defined. 

Proposition 2. If the closed path {a}x is fully retracing then V{{o\k) is a positive operator, else it's eigenvalues 
are either real or appear as complex conjugate pairs. 

Proof. This follows from two observations: 

S 2 n(p^)Si = tt(p£)*, (16a) 

K{pl{) = K(p$)\ (16b) 

Considering a matrix element of S^.1Z{p^i)Si: 

\S2'Tl{p2l)Sl) — (S2)ij;pq('R>(p2l))pq;rs(Sl)rs;kl = 0~iqdjp(TZ(p2i))pq;rs3rlO~sk 

V / ij-kl (17) 

= (R-(P2l))ji\lk = (P2l)jkM = {P*2l)il;jk = C^-(P2l))i r ,kh 



where the second last equality follows from the Hermiticity of p2i- This proves Eq. (|I6a[ ). using which and inserting 
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FIG. 2: (Color online) Eigenvalues of the matrix corresponding to P({1,2, 3}) for 100 random tripartite pure states with 
subsystem dimensions di = d,2 = d,3 = 10. The states are sampled according to the uniform Haar measure. The inset show an 
enlarged view of the region near the origin. 



Sf, i = 2, . . . , K as shown below, leads to 

- • • (k{ p t 2 {))* = v{{a} K y. 

In other words "P({o}k) is unitarily equivalent to its complex conjugate through the swap S\. Thus it immediately 
follows that V({o,}k) has a real characteristic polynomial and the eigenvalues appear as stated. 

That if the path is exactly retracing, ^({ajif) is a positive operator follows from Eq. (|16bl) which in turn follows 
from: 

(npti)) =(wiW = (P2i); a;V = (p 12 )^i/w = (^(^))* fl =(np T ii)f ■ m 

□ 

Finally, some consequences for tripartite higher dimensional pure states are now discussed. The closed path connect- 
ing all susbsystems of a pure state in one l oop maybe of special interest. For example for three qubits, tr("P({l, 2, 3})) 
over such a path is the Kempe invariant [171 . l24j . Its generalization to qudits then is evident from the formulation 
present above and is tr('P({l, . . . ,N})) valid for N subsystems of arbitrary dimensions. For example the eigenvalues 

of V({1, 2, 3}) = TZ(p^)TZ(pj2)T^{P2i) are shown in Fig. ((2|) for the case of random pure tripartite states of 3 systems 
each of dimension 10. It is observed that one real eigenvalue per realization is significantly larger in magnitude com- 
pared to others. This is a robust feature that is present in many operators, including elemental ones such as TZ(pW)- 
The origin of this is not hard to understand if one assumes a diagonally dominant density matrix, as typical ones 
indeed are. 

It is also not hard to see that when the state IV'123) is bi-separable, in particular IV'123) = ^12)® Ix3)i tne eigenvalue 
spectrum of ^({l^S}) consists of only one non-zero real value equal to tr^p^) 3 ]- F° r tlie path 1 —> 2 — s- 1 — >• 2, 
the same state gives for another invariant I' 6 = tr(7 , ({l,2, 1,2})) = tr[72.(p^2 )^-(P2i )P = [ tr (Pi)] 2 7 this also being the 
only non-vanishing eigenvalue of V({1, 2, 1, 2}). Further for tri-separable pure states the eigenvalues of 'P({1,2,3}) 
and 7 , ({1,2,1,2}) are all except one whose value is 1. 

In conclusion this work has pointed to a close connection between local invariants and the two operations that have 
hitherto been used to reveal entanglement, namely partial transpose and realignment. Properties of the resulting 
operators have been discussed and some simple consequences for qudit tripartite systems have been pointed out. It 
is hoped that this will lead to the construction of useful multipartite entanglement measures, and an understanding 
of invariants. 
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